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Comment on "Influence of Noise on Force Mea- 
surement" 

In a recent Letter [l| , Volpe et al. describe experiments 
on a colloidal particle near a wall in the presence of a 
gravitational field for which they study the influence of 
noise on the measurement of force. Their central result 
is a striking discrepancy between the forces derived from 
experimental drift measurements via their Eq. (1), and 
from the equilibrium distribution. From this discrepancy 
they infer the stochastic calculus realised in the system. 

We comment, however: (a) that Eq. (1) does not hold 
for space-dependent diffusion, and corrections should be 
introduced; and (b) that the "force" derived from the 
drift need not coincide with the "force" obtained from 
the equilibrium distribution. 

The problem of what should be the "correct" stochas- 
tic calculus was tackled in the early 1980s. The consensus 
was that, for a model in the form of a stochastic differ- 
ential equation (SDE), the calculus to be used, e.g. in a 
simulation, is part of the model itself. Correspondingly, 
starting from measured data, what we observe is a distri- 
bution function; but in the absence of further information 
and/or specific models, we cannot infer the underlying 
stochastic calculus Q . For a continuous physical system, 
with noise of (inevitably) finite bandwidth, we expect the 
Stratonovich calculus to apply Q. 

From the SDE Eq. (3) of 

dz = f{z)dt + g(z)dW = ^^dt + ^/2DZ{z)dW (1) 
we obtain the family of Fokker-Plank equations 


dP{z,t) 
dt 


ag{z)g'{z) + --/{z)^P{z,t) 

(2) 

where g'{z) = dg{z)/dz, and a is or 1/2 for the Ito 
or Stratonovich stochastic calculi respectively. In an ex- 
periment the diffusion (related to g'^{z)) and the drift 
(related to f{z) + ag{z)g'{z)) can be measured. To in- 
fer a, however, additional information (e.g. knowledge of 
f{z)) is needed. From ([2]), the drift velocity is 


Vd — —r — f(z) 
dt ' 


ag{z)g'{z) = 


F{z) ^ JDAz) 
"f{z) dz 


(3) 


This relation does not coincide with Eq. (1) of flj be- 
cause the nonlinearity of the diffusion coefficient enters 
the drift. Hence, it is impossible to derive the "force" 
F(z) from a measurement of the drift velocity, as in Eq. 
(1) of fl|, where it was assumed that F{z) = 'y{z)vd{z). 

In l] the force is also computed from Fe{z) — 
—dU{z)/dz where the potential U{z) = ~-kBT\x\{P{z)) 
is obtained from the the equilibrium distribution P{z). 
From (HI) we obtain 


U{z) 


f{z) + {a- l)giz)g'{z) 
g^z)/2 
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FIG. 1: (color online) Force computed from a simulation of 
^ using the Stratonovich calculus. Consistent with the ex- 
periment, we took f{z) = Be"'^^ + C with B = 770 pN, 
C = -5 fN, A; = (18 nm)"\ and D±{z) = D^z/{z + a) with 
z in nm, a — 700 nm, Doo ~ ksT /Qn-qR, 2R = 1.31 nm, 
T = 300 K, 77 = 8.5 X 10"=^ Pa s. 


l{z) 


1 dU{z) 
"f{z) dz 


= f{z) + l)g{z)g'{z) (4) 


dD^jz) 


which 


Eqs. dSI and (g]) differ by -g{z)g'{z) = 
is exactly the experimental discrepancy reported in [l[; 
this difference is independent of a, i.e. independent of the 
stochastic calculus used to describe the physical system. 

As a demonstration, we simulated ([T|) numerically for 
the Stratonovich calculus with the same definitions as [l[ , 
computing the average forces from the drift and equilib- 
rium distribution of the time sequence z{t). Fig.[T]shows 
that the force from the drift {vd"f{z)^ Eq. circles) dif- 
fers from that from the equilibrium distribution {Fe{z)^ 
Eq. ([U, full curve). When the drift result is corrected by 
the additional term —g{z)g'(z), however, we recover the 
equilibrium distribution result. Thus the discrepancy re- 
ported in P] has nothing to do with different stochastic 
calculi: it is simply a consequence of having two different 
definitions of "force" . Neither of them corresponds to the 
true microscopic force, and they coincide only where the 
diffusion coefficient happens to be constant. 

We are grateful to Ping Ao for alerting us to this prob- 
lem and to Mark Dykman for valuable discussions. 
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